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ABSTRACT 


The electric field acceleration of electrons out of a thermal plasma and 
the simultaneous Joule heating of the plasma are studied. Acceleration and 
heating timescales are derived and compared, and upper limits are obtained on 
the acceleration volume and the rate at which electrons can be accelerated. 
These upper limits, determined by the maximum magnetic field strength observed 
in flaring regions, place stringent restrictions upon the acceleration 
process. The role of the plasma resistivity in these processes is examined, 
and possible sources of anomalous resistivity are summarized. The 
implications of these results for the microwave and hard X-ray emission from 
solar flares are examined. 

The major conclusions are: 

(1) The simple electric field acceleration of electrons is found, in 

agreement with Spicer (1983), to be incapable of producing a large enough 

electron flux to explain the bulk of the observed hard X-ray emission from 

solar flares as nonthermal bremsstrahlung . For the bulk of the X-ray emission 
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to be nonthermal, at least 10 oppositely directed current channels are 
required, or an acceleration mechanism that does not result in a net current 
in the acceleration region is required. 

(2) If the bulk of the X-ray emission is thermal, a single current sheet can 
yield the required heating and acceleration timescales, and the required 
electron energies for the microwave emission. This is accomplished with an 
electric field that is much smaller than the Dreicer field (Ep/E ~ 10-50). 

(3) The rise time of the nonthermal emission is determined by the time needed 
to generate the required number of runaway electrons, rather than the time 



needed to accelerate the electrons to the required energies, which is 
generally a much shorter timescale. 

(4) The acceleration of enough electrons to produce a microwave flare 
requires the resupply of electrons to both the current sheet and the runaway 
region of velocity space. The electrons may be supplied by either the inflow 
of electrons that are already in the runaway region, or by the collisional 
scattering of electrons into the runaway region from the bulk plasma. For the 
latter case, the number of accelerated electrons is sensitive to the plasma 
resistivity. 

(5) To obtain the required heating timescale and electron energies, the 
resistivity in the current sheet must be much greater than the classical 
resistivity of a tplO cm , T e -10 K plasma. A plasma density -10 cm is 
required in the flaring region, or, if the density in the current sheet is 
less than 10 ^ cm - ^ , the resistivity in the sheet must be anomalous. Jj|- 
driven anomalous resistivity cannot be responsible for the higher resistivity 
unless the plasma electrons can first be heated to a temperature that is at 
least an order of magnitude greater than the ion temperature. 

(6) "Inertial resistivity", which was applied to flare heating by Duijveman, 
Hoyng, and Ionson (1981), is shown to always be negligible compared to the 
classical Coulomb resistivity. 

(7) Whenever the current is inductively determined by the magnetic field 
through Ampere's law, the energy that goes into Joule heating will always be 
greater than the energy that goes into accelerated electrons. 

Subject headings : particle acceleration — Sun: flares 


IV 



I. Introduction 


The impulsive phase of solar flares is characterized by the production of 
energetic charged particles and the rapid heating of thermal plasma. DC 
electric fields provide the simplest and most direct means of accelerating 
electrons out of a thermal plasma. The energy released in flares is generally 
understood to derive from the energy stored in magnetic fields and their 
associated currents. Therefore, most solar flare models result in the 
production of DC electric fields or, at least, electric fields that are not 
rapidly alternating (see Heyvaerts 1981 for a review). Since there is a 
macroscopic current associated with the electric field (J**aE), there will also 
be simultaneous Joule heating of the thermal plasma. 

Specific requirements for the number and energy of energetic electrons 
produced during a flare, and the timescales involved in accelerating them, are 
provided by microwave and hard X-ray observations of flares. The microwave 
emission from flares is understood to be gyrosynchrotron radiation from 
electrons with energies of 100 keV or greater. The hard X-ray emission (£25 
keV photons) can be interpreted as being either thick-target bremsstrahlung 
from nonthermal electrons (thin-target radiation may also contribute to the X- 
ray emission, but the process is less efficient), or thermal bremsstrahlung 
from hot, impulsively heated plasma. The relative contribution of thermal and 
nonthermal emission to the flare radiation is not well determined. 

In view of the constraints from flare microwave and hard. X-ray data, and 
the expectation that macroscopic electric fields and currents play an 
important role in solar flares, it is of interest to study the electric field 
acceleration of electrons and the Joule heating of flare plasma without 
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recourse to a specific model for the generation of the currents and electric 
fields. The basic physics of the electric field acceleration of “runaway" 
electrons is developed in Section II. Fundamental timescales for the 
acceleration of electrons out of a thermal plasma are derived and compared. 
Joule heating timescales are derived in Section III , and the coupling between 
Joule heating and electron acceleration is examined. 

The resistivity (thermal collision frequency) of a plasma plays an 
important role in determining both the rate at which electrons are accelerated 
and the Joule heating rate. Hence, the role of anomalous resistivity in these 
processes is examined in Section IV. The requirements for simultaneous 
electron acceleration and Joule heating in solar flares are studied in Section 

V. The major results of the paper are discussed and summarized In Section 

VI. 

Although the application in this paper is to solar flares, for which good 
microwave and X-ray spectral and timing data are available, it should be noted 
that most of the results presented here are directly applicable to any dynamic 
system of currents and magnetic fields where particle acceleration is 
apparent. In addition to the obvious application to stellar flares in 
general, runaway electron acceleration may be important in supernovae, close 
binary star systems, active galaxies, and other astrophysical systems where 
steep magnetic field gradients and, hence, significant DC electric fields can 
arise . 
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II. Electroa Acceleration 

The response of a thermal plasma to an imposed electric field, £, is to 
set up an electric current of current density env^-a^, where n is the 
number density of thermal electrons, e is the magnitude of the electron 
charge, is the mean drift velocity of the current-carrying electrons, and o 
is the conductivity of the plasma. A fraction of the current-carrying 
electrons, however, those, with a velocity greater than a critical velocity, 
v c , will be freely accelerated out of the thermal distribution (Dreicer 1959, 
1960). This can be seen from the single-particle equation of motion for an 
electron in the plasma: 

m Tt" * ” eE ~ mvv ( v ) » (1) 

where m is the electron mass and 


v(v) 


4 ine InA 

2 3 
m v 


( 2 ) 


InA * ln(12irn^) - 23.2 +'ln[( - q n - - )~ 1/2 ( — y — ) 3/2 I 

10 cin 10 K 


(2a) 


is the classical Coulomb collision frequency (Spitzer 1962, Knoepfel and Spong 
1979), (Collisions with Z>1 ions are neglected. See Knoepfel and Spong and 
Spitzer for a more precise determination of v and InA. Anomalous collision 
frequencies are discussed in Section IV.) The Debye length ■ v e /m e “ 

(kT/47rne 2 ) l/ 2 , where v e = (kT/m)^ is the electron thermal velocity and 
^“(Airne^/m)^ is the electron plasma frequency. 
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The bulk of the electrons are accelerated by the electric field until the 


steady state eE=mv^v e (equation 1) is obtained, where 


v 

e 


v(v ) 
e 


3m 


InA 
'e A 


10.00 


n. 

, n 9 ~ 3 

10 cm 


)(- 


10 K 




(3) 


(In using the thermal collision frequency, it is assumed that v^<<v e .) 

Since J^nev^, the resistivity of the plasma (n ■ 1/a) is found to be 
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Note that the value of the resistivity given by equations 3 and 4 is a factor 
of 2.6 greater than the Spitzer-Harm resistivity, which averages over the 
electron distribution, accounting for the skewing of the distribution in the 
driving electric field, in addition to electron-ion and electron-electron 
collisions • 

As can be seen from equations 1 and 2, the frictional drag on the 
electrons decreases with increasing particle velocity. Hence, electrons in 
the initial thermal distribution with a high enough velocity will not be 
confined to the bulk current, but will be freely accelerated out of the 
thermal distribution. The condition for electrons to be in this "runaway" 
regime is, from equation 1, 


eE > mvv(v) 
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(5) 


We see from equation 5 that a thermal electron of velocity v g will run away if 



the electric field strength is greater than 
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Ep is called the Dreicer field. For smaller field strengths (E<E^) 
electrons will run away that have velocities greater than 
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It is easy to see from equations 4 and 6 that the Dreicer field 
electrical resistivity are related through the expression 

E D 

— <* nev . 
n e 


Since n =» E/J =* E/nev^, the following relationships are obtained: 
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e 




Hence, when E=Eq, v d = ’ v c =v e anc * ctie c ^ erTna ^ electrons are in the 

regime. 

A rough estimate of the fraction of the electrons in the initial 
distribution that are in the runaway regime is 
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(Kaplan and Tsytovich 1973, Papadopoulos 1977). The timescale for the 
formation of the runaway tail is on the order of the collision time at the 
critical runaway velocity, or 

t t - < ii > 
e 

Although equation 10 estimates the initial number of electrons in the 
runaway regime, it does not give the total density or number of runaways 
produced, since, in addition to the acceleration of the initial thermal tail 
of electrons, new electrons are continously scattered into the runaway region 
by collisions. The kinetic theory of runaway production yields the following 
result for the runaway rate (Kruskal and Bernstein 1964 , Knoepfel and Spong 
1979): 


T * 0.35v e <I% 3/4 e*p[-2 V ^% - |<^) 2 ] S 
e e e 


-1 


( 12 ) 


This result is accurate for v»v . The rate at which runaway electrons are 

v. C 

collisionally produced is, then, N co n = "7 n ^j> where Vj is the volume of the 
current channel in which the electron acceleration is occuring (i.e., the 
volume of the acceleration region). 

Before obtaining final expressions for ^ c0 ]_i and the associated 
timescale, it is important to note that the volume Vj is not arbitrary, since 
the total allowable current in Vj is limited by the induction magnetic field 
associated with it. In other words, observations place an upper limit on the 
magnetic field strength that is present in active regions, and the magnetic 
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field associated with the total current I must not be stronger than this upper 


limit. This point was made earlier by Colgate, Audouze, and Fowler (1977) and 
by Colgate (1978), who argued that the hard x-ray emission from flares cannot 
be nonthermal bremsstrahlung from a directed beam of electrons, because the 
magnetic energy in the beam's self-generated field would exceed the magnetic 
energy available in the flare region by many orders of magnitude. It was 
realized by others, however, that the beam could induce a cospatial return 
current in the thermal plasma that would result in a negligible net current, 
thus avoiding the generation of a large induction field (Hoyng, Brown, and van 
Beek 1976; Knight and Sturrock 1977; Spicer and Sudan 1983). The current of 
interest here is driven by an associated electric field, however, as is 
required for the acceleration of runaway particles, and, therefore, there can 
be no cospatial return current (see also Hoyng 1977, Spicer 1983). Hence, in 
the acceleration region the induction magnetic field cannot be reduced by the 
presence of a return current. 

For a current sheet of width w, thickness <5r ( 6r«w) , and length L, with 
the current I flowing in the direction of L, the maximum value of the 
induction magnetic field is, from Ampere's law, Bj=(2ir/c)(I/w) . With I = Jw6r 
* nev^wSr, Bj a, (2re/c)(nv^6r) . Requiring Bj £ B, the maximum magnetic field 
strength in the acceleration region, yields <5r £ ( c/2re)(B/nv ( j) . Using 
equation 9 for v^, the thickness of the current sheet is found to be 

dr < 808( t ^- 5 -)( - 9 ? - 3 )~ 1 (-y— )~ 1/2 (^-) 2 cm. (13a) 
10 cm 10 K e 

Since Vj ” wLdr = A<$r, the acceleration volume is found to be 
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ATI r* p 2 

V < -r- — ( — ) ( — ) (sheet geometry). (13b) 

J — 2irne v v 
e e 

The sheet area. A, is limited by the observed scale of the acceleration 
region. Note that, if the current is flowing on the surface of a cylinder, as 
is likely to be approximately the case in a magnetic flux tube, then A=irRL, 
where R is the radius of the cylinder (Sr<<R). 

It is interesting to compare equation 13b with the corresponding result 
for a filled cylindrical geometry (uniform current density filling a circular 
cross section - a current filament). For this case Vj=irR L and 
Bj a, (2/c)(I/R) j “(2re/c)(nv c jR) . Hence, requiring Bj ( B gives R£( c/2ire) (B/nv^) 
and 


TR r 2 r 4 

V J ~ 1~ 2 } (filled cylindrical geometry). (14) 

4irn e 4 e e 


Since two spatial dimensions are constrained by the induction field, as 
compared to just one dimension for the sheet geometry, the constraint on Vj is 
more severe for a (filled) current filament. 

Using equations 12 and 13b, the runaway production rate is found to be 

7 - v - _i 

ft , i < 2.49x10 ABv ( — ) f(— — ) electrons s , (15) 

coll ~ e v v ’ 

e e 


where f(v c /v e ) < 1, with f(v c /v e ) - 1 when v c * 1.32v g : 


f(l£) - 4.66 (Ii) 11/4 exp [-2 1/ ^) -^<l£) 2 ]. 

e e e e 


(16) 
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The timescale for the production of N runaway electrons is t^ = N/N, or, from 
equation 15, 
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In addition to the constraint on Vj, the induction field limit puts a 
strict upper limit on N. The current associated with the runaway electrons is 
I^^eN. Since the induction field associated with this current is 
Bj=»(2Tr/c)(I run /w) , requiring Bj < B gives 


^Max " 1 ° 30( “^ )( T0§~G ) electrons s' 1 - (IS) 

10 cm 

Note that this result is independent of the temperature, density, and 
resistivity of the plasma. I run is a part of the total current, I, and will, 
always be less than (or equal to) I if I is determined by B through Ampere's 
law (see Spicer 1983 for further discussion of this point). Hence, this 
maximum value of N corresponds to the extreme case of I run =I. This result for 
N Max leads to a minimum value for the timescale required to generate N runaway 
electrons: 
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(19) 


For a cylindrical sheet or filament geometry, L is larger, and t2 in 

tlaX 1' 

smaller, by a factor of it (with R replacing w in equations 18 and 19). 
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In addition to generating the runaway electrons, it is also necessary to 
accelerate them up to the required energies. Neglecting the collision term in 
equation 1 , the time required to accelerate an electron from v c to a velocity 
v is, for a constant electric field. 


T 

a 


m(v-v ) E_ v , 

(_»)(! 

eE E v v ' e 

e e 


(^) 2 (— - l£)v _1 . 

v v v e 
e e e 


( 20 ) 


(If the final velocity is relativistic, v in equation 20 is replaced by yv, 
where Y is the usual Lorentz factor. Fewer runaways are generated if v c is 
relativistic, and none are generated if v c >c - Connor and Hastie 1975.) An 
important constraint on the acceleration time is the distance over which the 
acceleration can occur. For the simplest case of a constant electric field, 
the acceleration distance is 
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c a 
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( 21 ) 


2 

where We is the energy of the accelerated electron and W c =mv£/2. If the 

available distance over which the electrons can be accelerated is L, then 
requiring x < L yields an upper limit on t : 

d a 
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This requirement also leads to a lower limit on the electric field strength: 
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It is interesting that if we also require E<Eg, so that the entire particle 
distribution is not in the runaway regime, a lower limit on the collision 
frequency in the plasma is obtained (equations 6 and 23): 
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If the collision frequency is classical, this put3 a constraint oh the density 
and temperature in the current channel (equation 3). A time varying electric 
field is considered in Section III . 

In closing this section, it is worthwhile to compare the runaway tail 
formation, particle generation, and acceleration timescales, t T , t N , and r a . 
The time required to accelerate N electrons out of a thermal plasma to a 
velocity v is determined by the longest of these timescales. It can be seen 
from equations 11 and 20 that - [ ( v/v c )-l ] T^. therefore, f a >T T whenever 
v>2v c . Since v»v c for most astrophysical systems, the timescale is 
generally not important. Hence, the required time is determined by either T a 
or tjj. The timescale r a is restricted by an upper limit (equation 22), while 
a lower limit on tjj was found (equation 19). For the reference parameters 
used here in deriving the numerical coefficients, which are typical of a solar 
flare region, is seen to exceed T a « Hence, t^ is likely to be the dominant 
timescale for most astrophysical systems. 


III. Joule Heating 

The purpose of this section is to derive an estimate of the time required 
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to heat a current-carrying region through current dissipation, rather than to 
deveLop a detailed model for specific conditions. Hence, an estimate of the 
time required to heat a volume V to a given temperature will be obtained, 
without considering the details of the dissipation process, or heat loss 
mechanisms that may eventually balance the Joule heating. This timescale will 
be compared with the runaway timescales derived in Section II. Models for the 
heating of non-flaring coronal loops through current dissipation have been 
developed by a number of authors (Tucker 1973, Rosner et al. 1978, Hinata 
1980, Benford 1983). Joule heating in solar flares has been studied by Spicer 
(1981 a,b) and by Duijveman, Hoyng, and Ionson (1981). 

Since the energy dissipated by a current density J is J*E erg cm - ^ s”*, a 
Joule heating timescale Tj= nkT/(J‘E) can be defined. Using J=*E/n together 
with equations 4, 6, and 9 gives 


- nkT 
J = J*E 


&V 1 

o e 




(25) 


Note that tj is always longer than the runaway tail formation timescale 
(equation 11) as long as v c >v g , since Tj/t-j “ v c /v g . The timescale tj is 
shorter than the runaway acceleration timescale, x a (equation 20), whenever v 

> K v c /v e^ +l ^ v c* 

The timescale Tj is generally only a lower limit on the actual heating 
time, tj, since the volume of plasma to be heated, V, is usually much larger 
Chan the volume of the current channel, Vj. The global heating of the plasma, 
neglecting losses, is determined by the equation (cf. Duijveman et al. 1981) 
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fkilv - W! 0»J. 


(26) 


The heat deposited in the current channel may be conductively or convectively 
transported to the larger volume. It is implicitly assumed here that the heat 
is distributed throughout the larger volume on a timescale that is less than 
or on the order of the heating timescale. Otherwise the X-ray rise time will 
be determined by the heat transfer timescale rather than by the Joule heating 
timescale. 

It is apparent from equation 26 that in order to correct for the larger 
volume to be heated, Tj. should be divided by the factor £ = Vj/V, the fraction 
of the total volume occupied by the current channel. If the density in the 
current channel, n, is different from the average density in the volume V, n v , 
then Tj should also be multiplied by the factor ny/n, where Tj is evaluated 
for the physical parameters within the current channel. (The density in an 
equilibrium current sheet may be larger than that in the surrounding medium, 
for example, so that pressure balance is maintained.) Using equation 13 for 
Vj, the heating timescale is found to be 
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(27) 


Since £ is typically much less than unity, this timescale is generally much 
longer than Tj. 

In order to obtain a better feeling for the derived timescales, the 
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remainder of this section will be devoted to examining some simple models for 
the time development of the Joule heating and electron acceleration. The 
plasma density, magnetic field, and the volume V will be taken to remain 
constant. The Joule heating terra in equation 26 can be written as 
(J*E)Vj=*VjE 4 '/ n. If E and Vj both remain constant, this term is proportional 
co T 3/2, 3 i nce p varies as T ~^/2 ( C h e slow variation of InA with T will be 
neglected). On the ocher hand, if the cross sectional area of the current 
channel is determined by the ambient magnetic field so that Vj scales with 
temperature as T - ^^ (equation 13), the Joule heating term is independent of 
T. For this case the total current, I, remains constant, while I is not 
conserved if Vj remains constant. If J instead of E remains constant, the 
Joule heating term ( J *E)Vj=Vj is proportional co T ^/2. Vj necessarily 
remains constant in this case, and the total current is conserved. Using 
equation 26, the thermal evolution of the plasma can be obtained for these 
three cases by solving the equation 


dT 2 ~s 
dt = 3tj * 


(28) 


where T=T/Tq, Tq is the temperature of the plasma at time t=0, and tj is 
evaluated for T=Tq (equation 27). The solution to equation 28 is 


T * T 0 ti+ §<l-s)|-] 1/(1 " s) . (29) 

J 

For the case of E and Vj constant (s“*3/2), an increase in temperature by 
a factor of 10 requires t=2tj. The temperature increases rapidly, within the 
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time t=*3tj, until the heating is balanced by losses. Since E is constant, the 
runaway acceleration time t a “T a (equation 20). The runaway production time, 
t N , is given by equation 17 with T=Tq if tjj<<tj. The dominant effect of 
increasing the temperature on t N is to increase f(v c /v e ), decreasing t N . If 
t N (T Q ) ^ tj and i. tj, the runaway production time will be on the 

order of tj. If t^CTg) and t N (T Max ) are both longer than tj, the actual 
runaway production time will fall becween C^(T Max ) and ^(Tq). 

When E and I are held constant (s^O), the temperature increase is more 
gradual. A factor of 10 increase in temperature requires t=*13.5tj and the 
temperature increases linearly. Hence, c j(Tg) can significantly underestimate 
the actual heating time. 

When J and I remain constant (s=-3/2), the temperature increase is even 
more gradual. An order of magnitude increase in temperature requires 
t=189tj. This slow heating occurs because of the t“ 3/2 dependence of the 

classical resistivity so that, as the temperature increases, the heating rate 
decreases. Since E^nJ is not constant, the acceleration time is also 
modified. Using equation 29 with s=-3/2 for the time dependence of the 

temperature in E=»nJ, and solving equation 1 without the frictional drag term, 
gives 


t « \ t + i) 5/2 - i] . (30) 

a 5 J 3tj 

When T a «tj, C a ST a * If T a »tj, t a.~ T a /tj)^ 2 T a . is interesting that, 

unlike for the case of constant E, t^ increases with increasing T. 


The simple models presented here primarily indicate the care that must be 



exercised in applying Che derived timescales to astrophysical systems. They 
indicate how the actual physical times can vary from the estimated timescales 
under different physical conditions. Keeping these considerations in mind, 
however, they do indicate the conditions that are required for electric field 
acceleration and Joule heating to be applicable to a given astrophysical 
system. 

IV. Anomalous Resistivity 

The plasma resistivity can be much higher than the classical value 
(equations 3 and 4) if low-frequency electrostatic turbulence is present in 
the current channel. The current density J will induce the growth of such 
turbulence if its drift velocity (v^) exceeds a threshold velocity that Is 
greater than the plasma Ion sound speed, c g =(kT e /m i )^2 , The generation of 
anomalous resistivity by currents flowing parallel to the ambient magnetic 
field has been reviewed by Papadopoulos (1977). The generation of anomalous 
resistivity by currents flowing across the magnetic field has been reviewed by 
Papadopoulos (1980) and by Huba (1983). The growth of low-frequency 
electrostatic turbulence may also be stimulated by temperature or pressure 
gradients in the plasma (cf. Morrison and Ionson 1982), or by the presence of 
streaming, nonthermal electrons in the plasma (Papadopoulos 1977). The 
computed thresholds and wave levels for the J||-driven instabilities, which 
are the better studied of the current driven instabilities, are summarized in 
this section. 

For a current flowing along the ambient magnetic field, the instabilities 
that can lead to anomalous resistivity are the electrostatic ion cyclotron 
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(EIC), ion acoustic (IA.) , and Buneman instabilities. The Buneman instability 
is not of primary interest here, since its threshold velocity is given by v^i 
2v e and, therefore, the instability only occurs when the entire thermal 
electron distribution is in the runaway regime. Computer simulations indicate 
that heating causes the instability to saturate when v e =v^. If the current is 
driven by a constant electric field the instability oscillates, with both v^ 
and v g increasing in steps, keeping v^^e* 

The instability thresholds for the EIC and IA instabilities have been 
computed by Kindel and Kennel (1971) and Morrison and Ionson (1982) and 
references therein. These thresholds are a sensitive function of the ratio of 
the electron temperature to the ion temperature, T g /T^. The thresholds are 
plotted in Figure 1, based upon the results of Morrison and Ionson. These 
authors have accounted for the skewing of the electron distribution function 
in the driving electric field, and their (steady-state current) velocity 
thresholds are somewhat higher (by a factor of approximately 2) than those of 
Kindel and Kennel and others. (T g is the temperature of the unperturbed 
Maxwellian distribution function.) The effect of runaway electrons has not 
been included, which is likely to further increase the threshold velocities. 
For 0.1CT e /T^<8 the EIC instability has the lowest threshold (cf. Kindel and 
Kennel). Otherwise, the IA Instability is the first to go unstable. When 

T e = * T i* v thr^ EIC ^ =0,8v e ^ =35c s^* when V 8T i» v chr (EIC, IA) =8c g . For v thr =c g , 
T e must be more than an order of magnitude larger than T^ , with the IA mode 
being the unstable mode. 

Since the wave turbulence has the effect of increasing the collision 
frequency in the plasma, both the Joule heating and the runaways are 
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affected. If the effective collision frequency due to the wave turbulence, 
simply scales with particle velocity in the same way as the classical 
collision frequency, so that v)= y e ff ( v g /v) 3 , the results of the previous 

sections carry over to the case of a turbulent plasma by simply replacing v g 
with v g f f. Such a velocity dependence is expected for IA turbulence (Kaplan 
and Tsytovich 1973, Papadopoulos 1977). The turbulent collision frequency 
will be assumed here to always have the simple v J velocity dependence. 
Hence, in the turbulent plasma, v e ^ v e ff , n e f f“ 4ltv ef f f =(m/e)v e v ef f , 

v® ff -(mv eff v e 3 /eE) li/ 2 , etc. 

The approximate relationship between the anomalous collision frequency 
and the energy density in low-frequency IA turbulence, is 



W 


IA 


' nkT 


ID . 

e 


(31) 


The actual turbulence level that will be reached for a given instability is 
not well established, and depends upon the details of the system. Simulations 
of the IA Instability typically yield an effective collision frequency that is 
on the order of 


^ - 10 -2 iD = 10 6 (- 

eff e 


9 -3 

10 cm 


) l/2 (- 


10 s 


=r>-\ 


-1 


(32) 


Hence, an effective resistivity that is as much as a million times greater 
than the classical resistivity may be obtained if the drift velocity of the 
primary current exceeds the threshold for the IA instability. 

The anomalous resistivity that can result from the EIC instability is 


13 



less well established than the result for the IA instability. Quoted here is 
the result of Ionson (1976), which has been found to be consistent with a 
computer sumulation of the instability (Pritchett, Ashour-Abdalla, and Dawson 
1981): 


iff - °- 2 °i 3 2xl ° 4 <ToS-g><- 


10 s 




-1 


(33) 


where 0^ is the ion gyrofrequency. If both the EIC and the IA modes are 
present, the IA mode will generally provide the largest contribution to the 
effective resistivity. 

If the primary current is stable to the EIC and LA instabilities, it may 
be possible for a significant level of low-frequency IA turbulence to be 
generated by the runaway electrons. If the runaway electron distribution 
developes a positive slope so that it is unstable to the growth of Langmuir 
waves, and the energy density in these waves (normalized to nkT) reaches a 
level on the order of (v /v ) , where v r is the mean streaming speed of the 
runaway electrons, the waves can nonlinearly collapse to shorter wavelengths, 
generating ion density fluctuations in the process (the oscillating two-stream 
instability - OTS). The effective collision frequency resulting from the 
presence of these density fluctuations has been estimated by Papadopoulos and 
Coffey (1974) to be 


OTS 

\J 

eff 


n v Av 


(34) 


where Av r is the velocity spread in the runaway electron beam. For Av r /v r = 
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Hence , it may be 


1/3, v r /v e - 10, and n r =» 10 -3 n, for example, = 10~ 2 o> e . 

possible to generate an anomalous resistivity that is comparable to that 
estimated for the IA instability (equation 33). The generation of anomalous 
resistivity by streaming relativistic electrons has been studied by Scott et 
al. (1980), with application to the heating of gas in clusters of galaxies by 
cosmic ray electrons. 

An additional form of resistivity, "inertial resistivity" (Speiser 1970), 
has been applied to the heating of plasma in solar flares by Duijveman, Hoyng, 
and Ionson (1981). The argument for inertial resistivity is that if the time 
required for thermal electrons to cross the current sheet, ~dr/v e , where <5r is 
the thickness of the sheet (cf. equation 13), is less than a collision time, 
\T*, the electrons will be accelerated by the electric field while they are 
within the sheet and a net current density J IN = (ne 2 /m)( 5r/v e )E will flow. 
Since this current is therraalized in a few collision times, there is a 
contribution to the Joule heating rate with an effective collision frequency 
of v a ff “Vg/Sr. (This assumes, of course, that the crossing of the sheet by 
the thermal electrons is not suppressed by the ambient magnetic field.) Since 
<5r £ 1 km for solar flare conditions (see equation 13a and Section V.a.), the 
ratio v eff/ v e can be quite large. This result is misleading, however, since 

is derived on the basis of the current density J^, which is quite small 
when compared with the 3teady-state current density J=E/h e , where n e is the 
classical (Coulomb) resistivity. It is easy to see that the ratio of these 
current densities is Jjjj/J = v e <$r/v a , and the ratio of the corresponding Joule 
heating rates per unit volume is v e ^XN^ v e 3 3 ( v e ,5r / v e ) • Since the ratio of 
the volume containing J^ to the volume containing J is V IN /V J 3 W r ’ the 
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2 2 

ratio of Che total heating rates is v e ^XN V IN^ v e^ V J s v e^ r/,v e* T ^e orl -Si na l 
condition for inertial resistivity to be important was \> e <5r/v e < 1, however. 
Hence, heating due to inertial resistivity will never be significant compared 
to the heating by the primary current. The timescale for J^ N to be 
thermalized is the same as the timescale for the steady-state current density 
J to be established (a few collision times). Therefore, even if a steady- 
state current density (J) has not yet been achieved, inertial effects cannot 
contribute to the heating of the plasma until the larger current density J is 
established. 

V . Application to Solar Flares 

12 

A typical microwave flare requires that 10 - 10 J electrons be 

accelerated to energies on the order of 100 keV. It is interesting to compare 

this number of electrons with the number that can be supplied by the original 

distribution of thermal electrons in the current channel, N = n(v>v c )Vj, 

without considering the resupply of particles to the system. The maximum 

number of electrons that can be accelerated is obtained when v,=v »v and 

dee 

n(v>v c )=n, giving (equation 13b) 

S ( f-) ■ 8 - 1,<1029 electrons. (35) 

e 10 cm 10 K 

This result, which is independent of the density and the resistivity of the 
plasma, corresponds to accelerating all of the electrons in the current 
channel. A sheet area of 10 lo cm , if spread out in a single sheet parallel to 
the solar surface, corresponds to a 14"xi4" sheet, as seen from the earth. 
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The number of accelerated electrons decreases rapidly with an increase in the 
ratio v c / v e (equation 10). Note that a much smaller limit on N is obtained if 
the current is assumed to be confined to a cylindrical volume (equation 14), 
rather than to a sheet geometry. 

The maximum observed magnetic field strengths in flaring regions are on 
the order of 1000 G. It is interesting that under plausible, although 
somewhat extreme conditions it is possible for enough electrons to be 
accelerated from the original plasma in the channel to produce a microwave 
flare. In general, however, additional time is required for more electrons to 
be supplied from the thermal plasma. These particles must be obtained from 
plasma outside of Che acceleration region. The rate at which runaways are 
produced is then determined by the collisional rate, N cq11 (equation 15), or 
the rate at which new thermal electrons with v>v c flow into the current 
channel, whichever is greater. The current density J can supply electrons 
with v>v c to the acceleration region at the rate n(v>v ( ,)v ( jw6r. If plasma 
is flowing into the current sheet from the sides, as in a region where 
magnetic field merging is occuring, electrons in the runaway regime can be 


supplied at the rate N^* 2n 0 (v>v c )v IN A, where n Q is the density of the plasma 
outside of the current sheet and v-^ is the inflow speed of the plasma. For a 

steady-state magnetic neutral sheet, for example, v^ N is estimated to be 

l/„ 

~0 .01-0 . lv^, where v^= , B/(4irm^n 0 ) / 2 is the Alfven speed (Priest 1981). N cq ^ 

• • • 

is generally larger than N^j. The ratio of N c0 ]_]_ to N^ is 
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(36) 


v 
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where equation 10 has been used for n Q (v>v c ). For v c /v e ■ 5 and the other 

• • t • 

parameters as shown, for example, N co j^/N^ = 5*10 . Hence, is greater 

than N co j_j^ unless Vjjj is small or the thermal collision frequency in the 

channel is much greater than 10 s ^ . Irrespective of whichever process 

determines the runaway production rate, it cannot exceed the upper limit given 

by equation 18, however. 

• • 

VThen N RA does not exceed N co ^, an additional requirement for the runaway 

« 

rata to be determined by N co jli Is that thermal electrons be flowing into the 

• « • 

channel at a rate that is greater than “ Hence > either Nj ■ nv^wdr or 

= 2n Q v-j- N A must exceed N cq -q. Multiplying equation 36 by the factor exp[- 

(v c /v e ) I2\ gives the ratio It is easily seen that, unless there Is 

• • 

little or no inflow from the sides of the channel, does exceed N co n« The 

• • l 

ratio N c0 ^^/Nj is equal to simply yL/v^. The timescale Y 1 represents the 
time required for all of the thermal electrons entering the current sheet to 
be scattered into the runaway regime, while L/v^ is the time required for an 

average thermal electron in the current density J to traverse the length of 

• • • • 
the sheet. ^ co ll^j w ^ en l/ v d < I • It is interesting that when S coU >>N J’ 

the current is dominated by runaway electrons so that there is primarily 

• • 

particle acceleration and minimal Joule heating. N ^ can exceed Nj when the 

resistivity in the channel is high and v c /v e is small. As can be seen from 

• • • 

equations 13 and 18, however, if Bj=B. In this case N ^ will always 

♦ « • 

be less than or equal to Nj, since N co jj cannot exceed N^ ax * Hence, as has 
been emphasized by Spicer (1983), the primary current will always exceed the 
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runaway current if J is determined by B through Ampere's law and, therefore, 
the production of runaways will be maintained by Nj. 

It has been noted (cf. Norman and Smith 1978) that since the fraction of 
electrons in the runaway regime depends only upon v c / v e i and, therefore, v^/v e 
(equations 9 and 10), if v^ remains constant and there is a change in the 
plasma resistivity, there is no change in the number of accelerated 
electrons. Since the number of accelerated electrons depends in general upon 
the runaway production rate, N, rather than simply upon the fraction of 
particles initially In the runaway regime, however, this conclusion is not 
universally valid. If the runaway rate is determined by and v^ N is not 

sensitive to the resistivity in the channel, this conclusion remains valid. 
If the runaway rate is determined by N co ]^, however, the number of accelerated 
particles is sensitive to the resistivity in the current channel. 

The upper limit on N is particularly important for Interpreting the hard 

X-ray emission from flares. Interpreting the £.25 keV X-ray emission from a 

35 

flare to be thick-target, nonthermal bremsstrahlung requires a minimum of 10 
electrons s - ^ to be accelerated. Since observed magnetic field strengths are 
limited to B£1000 G, equation 18 shows that such a high particle flux cannot 
be achieved. Hence, the simple electric field acceleration of runaway 
electrons cannot produce a high enough electron flux to explain the bulk of 
the observed X-ray emission as non-thermal bremsstrahlung. This result agrees 
with a similar result obtained by Spicer (1983; see also Hoyng 1977). The 
hard X-ray emission could all be non-thermal if the flaring region is 
filamented into many small, oppositely directed (so that the net current is 
small or vanishes) current channels (with each channel conforming to the upper 
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limit on N). This requires _at_ least 10 individual current channels to 
produce a (i25 keV) hard X-ray burst, however. The plausibility of producing 
this situation, and its stability, requires further study. Alternatively, an 
efficient statistical acceleration process, which does not require or result 
in a net current in the acceleration region, may be able to generate the 
required flux of nonthermal electrons. Otherwise, although the highest energy 
X-rays may be nonthermal, the bulk of the hard X-ray emission must be thermal. 

In view of the preceding results, it is of interest to determine what 

conditions are required for the production of an impulsive hard X-ray burst 

through Joule heating, and an accompanying microwave burst through the 

acceleration of runaway electrons when a single current sheet is present. The 

microwave and hard X-ray rise times in solar flares range from seconds to 

hundreds of seconds (cf. Wiehl et al. 1983). For concreteness, a timescale of 

30 sec will be taken here as typical of both the microwave and hard X-ray rise 

times. The hard X-ray rise time is related to tj and, considering the 

discussion at the end of Section II, the microwave rise time is expected to be 

related to t^N/N. The exact relationship between tj and t^ and the rise 

times depends upon the individual flare, but for the purposes of this analysis 

it is reasonable to take tj=t^=30 sec. In equation 27 for tj, ny and V are 

not independent, since they are related through the (observed) emission 

2 

measure of the thermal X-rays, EM^n^ V. It is apparent from inspection of 
equation 27 that tj can only be of the right magnitude if v g is much greater 
than 10 s -1 in the current channel, or if the X-ray emission volume is much 
smaller than 10^ cm^ (and, therefore, is larger than 10^ cm - ^). Note also 
that B must be somewhat larger than 100 G, or w larger than 10 cm, if the 
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12 

lower limit on t^ (equation 19) is to be satisfied when N=10 electrons. 

Taking the runaway production rate to be determined by N co j_i> so that t N 
is given by equation 17, estimates for the required values of the collision 
frequency and v c /v e can be obtained. The value of v c /v e depends only upon the 
ratio of tj to tjj and the value of tiyV/N. For tj=t N and ny=10^ cm - ^ , V=*10 ^ 
cm^, and N^IO^, so that nyV/N =* 10^ (Case I), v c /v e =6 and, using equation 9, 
v^^l.ZCg. When n^lO^ cm - ^ and V=10^ cm^ , so that n^V/N ■ 10^ (Case II), 

v c =4v e and v^^.SCg. Similar results are obtained for v c if t^ is determined 

• 8 ~ — 1 * 
by with ^=30 sec and v IN =10 cm s 1 , rather than by N 

From equation 27 with tj a 30 sec, the required collision frequency in the 
current channel for Case I is found to be -10^ times greater than the 

Q 7 

classical collision frequency when n=10 cm and T=*10 K (this factor 
decreases somewhat if A or B is increased). Therefore, either the density in 
the current channel must be -10^ cm - ^ (n can be smaller if T is lower), or 
the resistivity in the channel must be anomalous. For Case II with iF'ny, the 
required collision frequency is only -10 times greater than the classical 
value. Densities as high as 10^ cm - ^ (but with T-10^ K) have been predicted 
for pre-flare current sheets (Syrovatskii 1976). Densities of 10 cm are 
commonly deduced from UV observations of flare loops (Dere et al. 1979 - it is 
not at present clear whether these densities are representative of the 
impulsive hard X-ray region, however). Hence, the impulsive microwave and 

hard X-ray emissions can be generated without the presence of anomalous 
resistivity within the current sheet. This is most easily accomplished with a 
plasma density -10^ cm ^ in both the current sheet and the X-ray emitting 
region. For a density of 10^ 4 cm~^ and T-10^ K, the thickness of the current 
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sheet is £1 an (equation 13a) and is less than or on the order of the 
gyroradius of a 100 keV electron (depending upon the value of B). It is also 
worth noting that for this case, with T-10^ K the bremsstrahlung cooling time 
in the sheet is on the order of the Joule heating time. If the density in the 
current channel is —10^-10*® cm - ^ , anomalous resistivity is required. For 
J| |-generated anomalous resistivity, Figure 1 and the derived values for v^ 
indicate that the ion acoustic instability will dominate. A high value of 
T e /T^ (> 10) is required for the instability, however. 

Although it is beyond the scope of this paper to consider heat transport 
mechanisms in flares, it is desirable to check the plausibility of the 
assumption made in Section III that the heat transfer timescale is less than 
or on the order of the Joule heating timescale. Taking the heat to be 
transported primarily in the direction perpendicular to the plane of the 
current sheet, the effective speed at which the heat must be transported is on 
the order of V/2Atj. This gives *-100 km s - ^ for Case I and -1,000 cm s - '’ for 
Case II. Both velocities are well below the Alfven speed, and the highest 
(Case I) is on the order of the ion sound speed. Hence, it Is reasonable to 
take tj to be the dominant timescale. It is certainly possible that in some 
flares the X-ray rise time is determined by the heat transport timescale, 
rather than by tj. In the absence of an additional heating mechanism, 
however, these flares will be characterized by a declining temperature during 
the rise phase of the flare, as the heat energy is spread over the larger 
volume. In this case, the results of the preceding paragraph are minimum 
requirements for the Joule heating to be accomplished with a single current 
sheet . 
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The time required to accelerate an electron from v c to an energy of 100 
keV (v =» 1.64x10*® cm s - *) is given by equation 20. For Case I T a ~10 - ^ sec, 
while for Case II t,-10 sec. These acceleration times are much less than 

a 

tfl, as expected. The energy that can be gained by an electron in a current 
channel of length L is found from equations 21, 9, and 6 to be 

W -W = 7.0 (— 5—: ) 1/2 ( S_-)(— )(-^)' 2 keV * (37 ) 

1 10 K 10 s 10 cm e 

For Case II with T=10^ K and v c / v e *4, W c =8 keV . Electron energies ~400 keV 

q 

can be attained in a distance of 10 cm for Case II. Much higher energies, 
~20 MeV , are possible for Case I. If the resistivity is due to a high density 
rather than to anomalous resistivity, however, the small thickness of the 
current sheet is likely to limit the electron energy that is attainable. It 
is interesting that, as was found for the Joule heating and runaway production 
timescales, the acceleration of electrons to an energy of 100 keV requires a 
collision frequency that is much greater than 10 s~*. 

VI. Discussion and Conclusions 

The rate at which the current I puts energy into Joule heating is ~TEL, 

while the rate at which energy goes into runaway electrons is I run E integrated 

over the length of the current sheet. Therefore, the energy that goes into 

particle acceleration will always be less than the energy that goes into 

• • 

heating if I run <l or» equivalently, if N<Nj. If I is determined Inductively 

• • • • 

with Bj=B, then, as was seen in Section V, Nj“N Max and N<Nj. Hence, the 
energy that goes into accelerated electrons will always be less than the 
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energy Chat goes into heating the plasma for such an inductively generated 
flare. This general result is in agreement with the conclusions reached by 
Smith (1980) from his analysis of specific inductive acceleration mechanisms. 

The runaway production rate chat is required for a microwave flare is 
Nh-(10^ electrons)/ (30 sec) * 3x10-^ electrons s - ^ - . Since this is on the 
order of N^ ax (equation 18), solar flares will have I run ~I when the electrons 
exit the acceleration region. The fact that the rate of electron acceleration 
required for a microwave flare is on the order of N^ ax may be an indication 
that an inductive flare model, with I determined by B through Ampere's law, is 
indeed correct. 

The direct acceleration of electrons by an electric field has the 
important consequence that the electron flux from one or a small number of 

current channels cannot be large enough to explain the bulk of the observed 

hard X-ray emission from flares as nonthermal bremsstrahlung . Hence, an 
important test of flare models involving the directed acceleration of 

electrons is to determine whether or not the lowest energy hard X-rays (20 to 
*•100 keV) from a given flare must be primarily nonthermal in origin. The 
current that is required to produce a typical hard X-ray burst has an 

induction magnetic field of at least 10^ G associated with it. Since such 
fields are not observed on the sun, if nonthermal electrons are required for 
these lowest energy X-rays, either the acceleration region must contain at 
least 10^ individual, oppositely directed current channels, or the 
acceleration mechanism must not require a directed current of electrons in the 
acceleration region. 

If the bulk of the hard X-ray emission is thermal, a single current sheet 
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is able Co generate Che required microwave and hard X-ray emissions. A 

9 9 

current sheet with dimensions —10 *10 cm that is inductively generated by a 
magnetic field of a few hundred Gauss can accelerate the required number of 
runaway electrons and deposit the required amount of heat on a timescale of a 
few tens of seconds if the density in the sheet and in the more extended X-ray 
emitting volume is -10^ cm - ^ . The thickness of the current sheet will be 
-400 cm, and the volume of gas to be heated is -200 times larger Chan the 
volume of the current channel. Electrons can be accelerated up to energies 

_ 9 

-400 keV in a timescale -10 sec. If the density in the current sheet is 
less than -10 cm , the resistivity in the sheet must be anomalous. A large 
collision frequency decreases the Joule heating timescale and the runaway 
production and acceleration timescales, and increases the maximum energy of 
the accelerated electrons. For an effective collision frequency -10^ s - * , 
electron energies up to -20 MeV are obtained. The required drift velocity of 
the current-carrying thermal electrons in the sheet is not high enough to 
expect J||-driven anomalous resistivity unless the thermal electrons are first 
heated to a temperature that is well above lOT^. If such an electron 
temperature can be attained, the ion acoustic instability will dominate. 

In a recent study of the hard X-ray and microwave emissions from a flare 
that occured on June 25, 1980, it was concluded that the emissions at the time 
of the peak 6 cm microwave radiation were most likely nonthermal (Holman, 
Kundu, and Dennis 1984). Since the power flux of electrons required for the 
hard X-ray emission (assumed to be thick-target bremsstrahlung) did not 
increase at this time, however, it was also concluded that either the 
acceleration mechanism must have only operated on already energetic particles 
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at this point in the flare, or two populations of energetic electrons were 
present. In view of the results obtained here, the best interpretation is 

that the bulk of the hard X-ray emission was thermal bremsstrahlung , while the 
secondary peak at the time of the microwave peak (the X-ray spectrum flattened 
at this time) corresponds to the acceleration of electrons to higher 
energies. The total electron power flux did not increase at this time 

because, although electrons were accelerated to high energies, the total 
energy in accelerated electrons was small compared to the energy in the hot 
thermal plasma. 

Some of the results presented in this paper have been summarized 
elsewhere (Holman 1984). These results will be applied to a specific flare in 
a subsequent publication (Holman and Kundu, in preparation; see also Holman 
1984). 
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